Modeling functional piezoelectricity in perovskite superlattices with competing 

instabilities 
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Based on the locality principle of insulating superlattices, we apply the method of Wu et al 
[Phys. Rev. Letter 101 , 087610 (2008)] to the piezoelectric strains of individual layers under fixed 
displacement field. For a superlattice of arbitrary stacking sequence an accurate model is acquired 
for predicting piezoelectricity. By applying the model in the superlattices where ferroelectric and 
antiferrodistortive modes are in competition, functional piezoelectricity can be achieved. A strong 
nonlinear effect is observed and can be further engineered in the PbTiOa/SrTiOa superlattice and 
an interface enhancement of piezoelectricity is found in the BaTiOs /CaTiOa superlattice. 

PACS numbers: 77.22.-d, 77.22. Ej, 77.80.-e, 77.84.Lf 
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Multicomponent ABO3 perovskite superlattices(SLs) 
provide a very promising way to design novel materi- 
als with multifunctional properties for device applica- 
tion Bridged by the interfaces, distinct instabilities 
belonging to individual bulk constituents are in strong 
competition with these artificial materials. Functional 
properties such as polarization |lU3] , piezoelectric Q , 
multiferroic 0^^: and dielectric responses are found to 
be highly sensitive to the interactive instabilities, some- 
times, resulting in unexpected enhanced functionalities. 

Paraelectric(PE)/ferroelectric(FE) SLs with both anti- 
ferrodistortive(AFD) and FE instabilities have attracted 
intense attention recently 0, 0, ■ The zone-boundary 
nonpolar AFD mode associated with the oxygen octa- 
hedral rotation and the zone-center polar FE mode are 
usually found to be exclusive to each other in bulk per- 
ovskite. However in SLs of PbTi03(PT)/SrTi03(ST) 
and BaTi03(PT)/CaTi03(ST), first-principles calcula- 
tions revealed that the AFD and FE can coexist with an 
interface reconstruction 0, Hi- In both cases, the com- 
peting AFD and FE modes at interfaces are predicted 
to enhance the polarization, which is consistent with the 
experimental observation [2, 11]. 

Clearly, for these insulating PE/FE SLs, both ionic 
and electronic interface effects should be well localized. 
The local electrostatic property will only be different 
from the bulk material within a few layers away from 
the interface under the fixed displacement field. Thus 
the interface dipoles can be accurately described by a 
truncated cluster expansion(TCE) model developed by 
Wu et al. 12| ; in which the electronic states of interfaces 
are represented by the maximally localized Wannier func- 
tions th rough a unitary transformation from Bloch-like 
orbitals jl3l-[l5|. Based on the above, SL design can be 
performed in the AFD/FE competing system, where the 
interface is expected to increase the overall FE. For quite 
some time, the focus was to explore the functional po- 
larization ji-HITel-flill. Piezoelectricity [Til-E^. another 
important functionality, describes the coupling between 
polarization and strain. Its interface mechanism and SL 



design rule have not yet been addressed in these intrigu- 
ing systems with competing instabilities. 

In this letter, we show that functional piezoelectricity 
can be designed in the SLs through the AFD and FE 
competition and its electric field dependence. We fur- 
ther apply our TCE model to the piezoelectric strain of 
the SLs in a fixed displacement field. In combination 
with a similar treatment of Wannier based layer polar- 
ization, we arrive at an accurate modeling for predict- 
ing the piezoelectric tensor for an arbitrary sequence of 
SL. In bulk ST, we discover a strong nonlinear piezoelec- 
tric response originating from the completely suppressed 
AFD in a large applied electric field. We use the model to 
demonstrate that the electric field of this nonlinear piezo- 
electricity can be reduced to a much smaller magnitude 
with an largely increased field tunability by changing the 
PT fraction in PT/ST SLs. With the model, we are able 
to systematically study the interface effect on the piezo- 
electric response in both BT/CT and PT/ST SLs. In 
BT/CT SLs, we find a novel interface enhancement of 
piezoelectricity. 

The first-principles calculations are carried out in two 
series of SLs nPT/mST and nBT/mCT both with AFD 
and FE instabilities. These include the bulk PT, ST, 
BT, CT and all period-four SLs(n + m = i) stacked in 
the [001] direction. The fixed in-plane lattice constant 
is chosen as 00=7.275 bohr, the computed equilibrium 
lattice constant for bulk cubic SrTi03. We assume the 
epitaxial growth of the SLs on SrTi03 which is consistent 
with the recent experiments 0, Q . Both the FE mode 
along [001] and the AFD mode associated with oxygen 
octahedral TiOe rotation around [001] axis have been 
explicitly included in the sub-lattice structure of PAbm 
space-group symmetry. To this end, a 40-atom tetragonal 
supercell is adopted with lattice vectors of length V2ao 
in the [110] and [110] directions and [110] directions and 
c ~ 4ao along the [001] direction. 

We use density-functional theory implemented in the 
LAUTREC code package Q to perform structural relax- 
ations and electron minimizations at fixed electric dis- 
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FIG. 1: (color online) Dependence of Ti02 (T) (a) layer height 
(b) layer polarization(relative to the average of the TiOa 
planes of bulk PbTiOa and SrTiOa) on chemical environment 
in PT/ST SLs. 



placement fields. Local-density approximation ^24] is 
adopted and plane-wave calculations are implemented 
in the projector augmented-wave framework [25| . We 
used a plane- wave cutoff energy of 80 Ry and a 4x4x1 
Monkhorst-Pack k mesh. 

Our model for piezoelectricity starts from the decom- 
position of the SLs into AO and BO2 layers along [001] di- 
rection. For each fixed-D field, both Wannier-based layer 
polarization Pi{D) fl^] and the layer height hi{D) Q 
(which is directly related to piezoelectric strain rji{D)) of 
each individual layer are computed from the relaxed elec- 
tronic and ionic structure respectively. Working with the 
constrained-!? field framework 23|, we employ longitu- 
dinal boundary condition which limits the force constant 
matrix to short range interaction of a few neighbors. 



In Fig. 1, we present the representative hj{D) and 
Pj{D) curves in PT/ST SLs. As expected, the local- 
ity principle is satisfied on not only the layer polariza- 
tion [l3| but layer height. As a result both hj{D) and 
Pj{D) curves separate into clusters depending on the 
nearest neighbor chemical environment. In can be seen 
that Pj{D) is mostly determined by the identity of near- 
est neighbors and has a much weaker dependence on the 
second neighbors. Compared with Pj{D), hj{D) shows 
a similar behavior with an even stronger localization. A 
clear local inversion symmetry breaking is also present 
in hj{D) as well as pj{D). As an example, the inversion 
symmetry breaking introduced by the first neighbors has 



a large asymmetric behavior for the hj{D) of Ti02 in the 
middle of the S(T)P sequence. This asymmetric behavior 
becomes much weaker when the local inversion symmetry 
breaking occurs only on the second neighboring layers in 
SS(T)SP. In previous work [l^, it has been shown that a 
TCE model accurately captures the dependence of pi (D) 
on its local compositional environments. The similar lo- 
cality of hi{D) indicates that hj{D){ds well aspj{D)) can 
be accurately described by the TSE as 

His}) = Jo + Y,iJh^s, + J[^^s^^) 

■i 



ijk 



(1) 



Our TCE model will include the cluster interaction of 
one particular layer with the neighboring layers up-to sec- 
ond nearest neighbors and only include up to two-body 
interaction terms. As a result, all the truncated clus- 
ter terms are effectively included in the compositional 
environment of the period-four SLs. All the computed 
Pj{D) and hj{D) will serve as a first-principles database. 
The J terms are D-dependent effective cluster interac- 
tion coefhcients and are computed from the database. 
The "pseudo-spin" variable is defined as Si = 1,0 and is 
used to identity the «th AO layer as either PbO(BaO) or 
SrO(CaO) respectively. The expressions for /iao, ^Ti02i 
PAO and PTi02 with all the fitted ECIs will be given 
in the on-line supporting material [26| . The total su- 
percell lattice constant is obtained from the sum of all 
the individual layer heights as h{D) = ^ihi{D) and 
the piezoelectric strain of the SL can be calculated by 
r]{D) = {h{D) - h{D = 0))/h{D = 0). The electric equa- 
tion of state of ry as a function of electric field f (2?) will 
be obtained by numerical inversion The electric field is 
computed by £{P) — D{P) — 4:TtP, where the total polar- 
ization of the SL is given by P{D) = h{D)-^ Ej ft(^)- 
It is then straightforward to compute the piezoelectric 
coefficient ^33 = drj^s/dS. 

Here we want to stress that the modeling of the indi- 
vidual layer hight hj{D), instead of the total height h{D) 
of the SL, is important for a more accurate piezoelectric 
model. This is because the piezoelectricity describes the 
derivative of the strain with respect to the electric field, 
which is sensitive to the accuracy of the strain model. 
Due to the short-range nature, the models of Pj{D) and 
hj{D) can be used to compute r]{D) and f (2?) for an 
arbitrary sequence of SL. This enables us to study the 
piezoelectricity resulting from the AFD/FE competition 
under an applied electric field or at the interfaces sys- 
tematically. 

It is well established that AFD and FE will compete 
with each other in bulk ABO3 [13, IHliil. Under an 



applied electric field, the FE will be strengthened by the 
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FIG. 2: (Color online) TiOe rotation as a function of D field 
for TiOe octahedral in bulk ST, strained bulk PT and strained 
bulk CT(shown in insert). 

coupling between the cleetric field and polar mode and 
the SL will be further polarized along the field direction. 
It is thus expected that AFD will be suppressed with 
increased field. Indeed in Fig. 2, we can see that the 
octahedral rotation is largely suppressed with increased 
D magnitude in both bulk ST and strained bulk PT. 
Surprisingly at D ^ 0.4 C/m'^ and D ^ 0.6 C/m?, the 
octahedral rotation will be completely suppressed result- 
ing in a FE only phase. It is also important to notice 
that the TiOg rotation in bulk PT is hidden and can not 
be observed experimentally. This is because the AFD 
instability is completely suppressed in the metastable 
region before the spontaneous polarization is reached 
(D ~ P = 0.8 C/to2). In the PT/ST SLs, the spon- 
taneous polarization will be reduced and the AFD rota- 
tion will be recovered in the PT fraction. In contrast, 
octahedral rotations in strained bulk CT have a much 
weaker field dependence. The TiOg rotation in CT is 
only slightly decreased from 13° at Z? = as shown in 
the inset of Fig. 2 This is consistent with the much 
stronger AFD than FE instability which results in the PE 
ground state for bulk CT. On the other hand BT shows 
no AFD rotation at all for the whole field range. This is 
expected from the fact that bulk BT highly resists AFD 
rotation with a robust FE ground state. 

Since ST is PE at ground state, the disappearance of 
the AFD instability indicates a physical phase transition 
driven by an applied electric field. During the phase tran- 
sition, SL structure changes from AFD/FE into a FE 
only phase accompanied by a structure softening, which 
is signaled by a strong nonlinear piezoelectric at finite 
electric field. One can clearly see an additional peak of 
d^ziS) in bulk ST in Fig. 3. However, it occurs at a very 
high electric field(centered around 600 Mv/m). This is 
probably the reason why this phase transition has not 
been addressed yet in the literature. 

Under the fixed-D field, both the FE and AFD or- 
dering is short-ranged. So it should be understood that 
the D field is the fundamental variable that drives the 
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FIG. 3: (Color online) Dependence of d33(£) piezoelectric 
coefficient on PTO fraction in the nPT/mST SLs. 



phase transition at which the threshold D field is reached. 
Keeping the above concept in mind, we propose that D- 
field can be increased by a highly polarizable PT com- 
ponent in the PT/ST SLs. As a result of the modified 
energetics of AFD and FE, the above phase transition, 
as well as the piezoelectric anomaly, can be observed at 
a much lower electric field. We then use our developed 
model to test the above idea. We set up a nPT/mST 
SL where n + m = 200 and gradually increase the PT 
fraction from ST bulk n — 0. The model prediction of 
the piezoelectric coefficient is presented in Fig. 3. As ex- 
pected, we can see that the SLs become more polarized 
as the PT fraction is increased. The SL starts from a 
PE state and becomes FE after the PT fraction in the 
SL is larger than ~ 30%. As a result of increased FE or- 
dering, the piezoelectric coefficient is generally increased 
with the anomalous increase corresponding to the phase 
transition in the ST component. Furthermore, the center 
of this anomaly is gradually shifted to lower electric field 
(centered at 55 mV/m for 60% PT fraction). In having 
increased polarizability from a higher PT fraction in the 
SLs, the width of the piezoelectric anomaly also becomes 
smaller indicating the greatly increased tunability of ap- 
plied electric field. 

AFD and FE competition at interfaces is also found 
to be of critical importance for its functional proper- 
ties [1, 01 • Here, focusing on the piezoelectricity, we 
can use our model to study the interface effect in both 
BT/CT and PT/ST systems. To the above end, we set 
up the SL model of nPT/nST and nBT/nCT and grad- 
ually increase n from dense interface limit (n = 1) to the 
interface free limit {n = oo). The model prediction of the 
d^zi^) is presented in Fig. 4 (the predicted spontaneous 
polarizations Ps as a function of n are also plotted in the 
inserts of Fig. 4(a) and (b)). The d33(f ) from directed 
first-principles are also given for comparison. One can see 
that our model is very accurate in reproducing the first- 
principles results. In both systems the SL will have the 
largest at the dense interface limit {Ps= 0.187 C/m^ 
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FIG. 4: (Color online) Dependence of d33(£) piezoelectric 
coefficient on the density of interfaces in (a) nPT/nST SLs 
(b) nBT/nCT SLs. 

for IBTICT and 0.327 C/m^ for IPTIST). The will 
monotonically decrease until saturate at 0.21 C/w? and 
0.14 C/tc? for nPT/nST and TiBT/nCT in the interface- 
free limit. This is consistent with the interface enhance- 
ment of polarization found in these materials. Surpris- 
ingly, we see also an interface enhanced piezoelectricity 
in nBT/nCT where c?33(f) gains its maximum magni- 
tude at n = 1. In contrast the nPT/nST SLs show the 
opposite effect, where the nonlinear effect corresponding 
to the ST phase transition starts to be observed when 
n = 3. 

The increased dj,^{£) indicates the structural soften- 
ing introduced by the interfaces. In FE materials, the 
structural softening can be found when the system is ap- 
proaching the FE/PE phase boundary e.g. with chang- 
ing in-plane epitaxial strain [7|]. As a result, the dielec- 
tric and piezoelectric responses diverge in the vicinity of 
phase boundary. Indeed, we do see a similar interface 
dependence of the static dielectric response (including 
the piezoelectric mediated component) ;26] confirming 
the additional softening of the polar mode due to the in- 
terfaces. Strikingly, in BT/CT SLs the polarization and 
piezoelectric response can be both increased by the in- 
terface which is crucial for SLs design for multifunctional 
properties. 

In summary, we have developed an accurate model that 
can the predict of piezoelectric coefficient for an arbitrary 
sequence of SL using the first-principles results of short- 



period SLs only. The predictive power of the model has 
been demonstrated in the PT/ST and BT/CT systems 
in which AFD and FE are in strong competition. Func- 
tional piezoelectricity can be designed in the PT/ST and 
BT/CT SLs under an applied electric field or a function 
of interface density. 
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